Leak-rate of seals: effective medium theory and comparison with experiment 
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Seals are extremely useful devices to prevent fluid leakage. We present an effective medium theory 
of the leak-rate of rubber seals, which is based on a recently developed contact mechanics theory. 
We compare the theory with experimental results for seals consisting of silicon rubber in contact 
with sandpaper and sand-blasted PMMA surfaces. 



o 
o 

(N 
> 

o 



O 

C3 



I 

-a 
c 

o 
o 



> 

On 

o 



ON 

o 



X 



1. Introduction 

A seal is a device for closing a gap or making a joint 
fluid tight [1]. Seals play a crucial role in many modern 
engineering devices, and the failure of seals may result 
in catastrophic events, such as the Challenger disaster. 
In spite of its apparent simplicity, it is not easy to pre- 
dict the leak-rate and (for dynamic seals) the friction 
forces [H- The main problem is the influence of surface 
roughness on the contact mechanics at the seal-substrate 
interface. Most surfaces of engineering interest have sur- 
face roughness on a wide range of length scales^, e.g, 
from cm to nm, which will influence the leak rate and 
friction of seals, and accounting for the whole range of 
surface roughness is impossible using standard numerical 
methods, such as the Finite Element Method. 

We have recently presented experimental results for 
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and compared the results 
I, H, Q , which is based on 
percolation theory_^and a recently developed contact me- 



13 . Il4j . Here we will 



the leak-rate of rubber seals p[ 
to a "single-junction" theory 

eoryand a recently c 
chanics theory 0,1, ©, M, [H, El 
report on new experimental data, and compare the exper- 
imental results with the single-junction theory, and also 
to an extension of this theory presented below, which is 
based on the effective medium approach. 

2. Single-junction theory 

We first briefly review the leak-rate model developed in 
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FIG. 1: Rubber seal (schematic). The liquid on the left-hand- 
side is under the hydrostatic pressure P a and the liquid to the 
right under the pressure Pb (usually, Pb is the atmospheric 
pressure). The pressure difference AP = P a — Pb results in 
liquid flow at the interface between the rubber seal and the 
rough substrate surface. The volume of liquid flow per unit 
time is denoted by Q, and depends on the squeezing pressure 
Po acting on the rubber seal. 
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FIG. 2: The contact region at different magnifications £ = 3, 
9, 12 and 648, is shown in (a)-(d) respectively. When the 
magnification increases from 9 to 12 the non-contact region 
percolate. At the lowest magnification f = 1: A(l) = Ao- The 
figure is the result of Molecular Dynamics simulations of the 
contact between elastic solids with randomly rough surfaces, 
see Ref. @|. 



Ref . [1, [|| . Consider the fluid leakage through a rub- 
ber seal, from a high fluid pressure P a region, to a low 
fluid pressure Pb region, as in Fig. [1] Assume that the 
nominal contact region between the rubber and the hard 
countersurface is rectangular with L v , with 

L y > L x . We assume that the high pressure fluid region 
is for x < and the low pressure region for x > L x . 
We "divide" the contact region into squares with the 
side L. x = L and the area A® = L? (this assumes that 
= Ly/ L x is an integer, but this restriction does not 
affect the final result). Now, let us study the contact 
between the two solids within one of the squares as we 
change the magnification £. We define £ = L/X, where 
A is the resolution. We study how the apparent contact 
area (projected on the xy-plane), A((), between the two 



solids depends on the magnification £. At the lowest mag- 
nification we cannot observe any surface roughness, and 
the contact between the solids appears to be complete 
i.e., A(l) = Aq. As we increase the magnification we will 
observe some interfacial roughness, and the (apparent) 
contact area will decrease. At high enough magnifica- 
tion, say £ = Cc, a percolating path of non-contact area 
will be observed for the first time, see Fig. [2] We denote 
the most narrow constriction along this percolation path 
as the critical constriction. The critical constriction will 
have the lateral size A c = L/Q c and the surface separa- 
tion at this point is denoted by u c . We can calculate u c 
using a recently developed contact mechanics theory [l2j 
(see below). As we continue to increase the magnifica- 
tion we will find more percolating channels between the 
surfaces, but these will have more narrow constrictions 
than the first channel which appears at £ = £ c , and as a 
first approximation one may neglect the contribution to 
the leak-rate from these channels^. 

A first rough estimate of the leak-rate is obtained by 
assuming that all the leakage occurs through the critical 
percolation channel, and that the whole pressure drop 
AP = P&—Pb (where P a and Pb is the pressure to the left 
and right of the seal) occurs over the critical constriction 
(of width and length A c s» £/£ c and height u c ). We will 
refer to this theory as the "single-junction" theory. If 
we approximate the critical constriction as a pore with 
rectangular cross section (width and length A c and height 
u c « A c ), and if assume an incompressible Newtonian 
fluid, the volume-flow per unit time through the critical 
constriction will be given by (Poiseuille flow) 



q = *1 AP, 



(1) 



where n is the fluid viscosity. In deriving (1) we have 
assumed laminar flow and that u c << A c , which is al- 
ways satisfied in practice. We have also assumed no-slip 
boundary condition on the solid walls. This assumption 
is not always satisfied at the micro or nano-scale, but is 
likely to be a very good approximation in the present case 
owing to surface roughness which occurs at length-scales 
shorter than the size of the critical constriction. Finally, 
since there are N = L y /L x square areas in the rubber- 
countcrsurface (apparent) contact area, we get the total 
leak-rate 



Q = 



L x 12 V 



AP. 



(2) 



Note that a given percolation channel could have several 
narrow (critical or nearly critical) constrictions of nearly 
the same dimension which would reduce the flow along 
the channel. But in this case one would also expect more 
channels from the high to the low fluid pressure side of 
the junction, which would tend to increase the leak rate. 
These two effects will, at least in the simplest picture 



where one assumes that the distance between the critical 
junctions along a percolation path (in the x-direction) is 
the same as the distance between the percolation chan- 
nels (in the y-direction) , compensate each other (see Ref. 
0). The effective medium theory presented below in- 
cludes (in an approximate way) all the flow channels. 

To complete the theory we must calculate the separa- 
tion u c of the surfaces at the critical constriction. Wc 
first determine the critical magnification ( c by assuming 
that the apparent relative contact area at this point is 
given by site percolation theory. Thus, the relative con- 
tact area A(()/Aq « 1 — p c , where p c is the so called site 
percolation threshold [20]. For an infinite-sized systems 
p c ss 0.696 for a hexagonal lattice and 0.593 for a square 
lattice 20]. For finite sized systems the percolation will, 
on the average, occur for (slightly) smaller values of p, 
and fluctuations in the percolation threshold will occur 
between different realizations of the same physical sys- 
tem. We take p c w 0.6 so that A(( c )/A w 0.4 will 
determine the critical magnification £ = £ c . 

The (apparent) relative contact area A(()/Aq at the 
magnification £ can be obtained using the contact me- 
chanics formalism developed elsewhere 0, H, E3, EL EH, 
where the system is studied at different magnifications £. 
We have@,H 
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where the surface roughness power spectrum 
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where (...) stands for ensemble average. Here E and v 
are the Young's elastic modulus and the Poisson ratio of 
the rubber. The height profile h(x) of the rough surface 
can be measured routinely today on all relevant length 
scales using optical and stylus experiments. 

We define ui(C) to be the (average) height separating 
the surfaces which appear to come into contact when the 
magnification decreases from £ to ( — A£, where A( is a 
small (infinitesimal) change in the magnification. In Fig. 
|3{a) the black area is the asperity contact regions at the 
magnification £. The green area is the additional contact 
area observed when the magnification is reduced to C~^C 
(where A£ is small) [2l[ . The average separation between 
the solid walls in the green surface area is given by ui(C)- 
Fig. El^b) shows the separation between the solid walls 
along the dashed line in Fig. [3ja). Since the surfaces 
of the solids are everywhere rough the actual separation 
between the solid walls in the green area will fluctuate 





FIG. 3: (a) The black area is the asperity contact regions at 
the magnification f. The green area is the additional contact 
area observed when the magnification is reduced to f — A( 
(where A£ is small) . The average separation between the solid 
walls in the green surface area is denoted by «i(C). (b) The 
separation between the solid walls along the blue dashed line 
in (a). Since the surfaces of the solids are everywhere rough 
the actual separation between the solid walls in the green area 
will fluctuate around the average Wi(C)- At the most narrow 
constriction the surface separation is u c . 
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FIG. 4: An asperity contact region observed at the magnifica- 
tion f . It appears that complete contact occur in the asperity 
contact region, but when the magnification is increasing to the 
highest (atomic scale) magnification £1 , it is observed that the 
solids are actually separated by the average distance u(Q. 
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FIG. 5: Effective medium theories take into account random 
disorder in a physical system, e.g., in a granular metallic film. 
The equation determining the "effective medium" (e.g., the 
effective conductivity cr e g) is obtained by calculating some 
properties of the effective medium and demanding that the 
same result is obtained by embedding in the effective medium 
a circular region of one component of the original system, and 
then averaging over the different component, with weights 
determined by the fractional areas of the various components 
in the original physical system. 



around the average iii(0- Thus we expect u c = crai(£ c ), 
where a < 1 (but of order unity, see Fig. [3jb))[22j]. We 
note that a is due to the surface roughness which occur 
at length scales shorter than A c , and it may be possible 
to calculate (or estimate) a from the surface roughness 
power spectrum, but no such theory has been developed 
so far and here we treat a as a fit parameter. 

ui(0 is a monotonically decreasing function of £, and 
can be calculated from the ave rag e interfacial separation 
w(C) and A(C) using (see Ref. 

Ul (0=u(0+u\C)A(C)/A'(C). 

The quantity u(Q is the average separation between the 
surfaces in the apparent contact regions observed at the 
magnification £, see Fig. 0] It can be calculated from(l^| 

dqq 2 C(q)w{q) / dp' _ e -K«.C)p /B ] ) 



6(0 = A 
where p(Q = P A /A(() and 
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The function P{q 1 p, Q) is given by 
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where s(q, Q = w(q, 0/23*. 
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3. Effective medium theory 

The single-junction theory presented above assumes 
that the leak-rate is determined by the resistance towards 
fluid flow through the critical constriction. In reality 
there will be many flow channels at the interface. Here 
we will use the 2D Bruggeman effective medium theory 
to calculate (approximately) the leak-resistance resulting 
from the network of flow channels. 
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FIG. 6: Experimental set-up for measuring the leak-rate 
of seals. A glass (or PMMA) cylinder with a rubber ring 
attached to one end is squeezed against a hard substrate 
with well-defined surface roughness. The cylinder is filled 
with water, and the leak-rate of the water at the rubber- 
countersurface is detected by the change in the height of the 
water in the cylinder. 



Using the 2D Bruggeman effective medium theory we 
get (see Ref. [Hj|, Fig. and Appendix A): 



Q = ^a cS AP, 



(3) 



where AP = P a — Pb is the pressure drop and where (see 
Appendix A) 

— =[ da P(a) 



A Q J er ff + o-(C) ' 



where 
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(4) 



(5) 



Eq. (4) is easy to solve by iteration. 

It is not clear that the effective medium theory is better 
than the single-junction theory. One problem with this 
theory is the following: In the effective medium model 
there is no correlation between the size of a region and 
the (average) separation between the surfaces in the re- 
gion. In reality, the regions where the surface separation 
is large form large compact (or connected) regions (since 
they are observed already at low magnification). 

4. Experimental 

We have performed a very simple experiment to test 
the theory presented above. In Fig. [5] we show our set- 
up for measuring the leak-rate of seals. A glass (or 



PMMA) cylinder with a rubber ring (with rectangular 
cross-section) attached to one end is squeezed against a 
hard substrate with well-defined surface roughness. The 
cylinder is filled with water, and the leak-rate of the fluid 
at the rubber-countersurface is detected by the change 
in the height of the fluid in the cylinder. In this case 
the pressure difference AP = P a — Pb = pgH, where 
g is the gravitation constant, p the fluid density and 
H the height of the fluid column. With ff a 1 m 
we get typically AP 0.01 MPa. With the diameter 
of the glass cylinder of order a few cm, the condition 
Po >> AP (which is necessary in order to be able to 
neglect the influence on the contact mechanics from the 
fluid pressure at the rubber-countersurface) is satisfied 
already for loads (at the upper surface of the cylinder) 
of order kg. In our study we use a rubber ring with 
the Young's elastic modulus E = 2.3 MPa, and with 
the inner and outer diameter 3 cm and 4 cm, respec- 
tively, and the height 0.5 cm. The rubber ring was 
made from a silicon elastomer (PDMS) prepared using 
a two-component kit (Sylgard 184) purchased from Dow 
Corning (Midland, MI) . The kit consist of a base (vinyl- 
terminated polydimethylsiloxane) and a curing agent 
(methylhydrosiloxane-dimethylsiloxane copolymer) with 
a suitable catalyst. From these two components we pre- 
pared a mixture 10:1 (base/cross linker) in weight. The 
mixture was degassed to remove the trapped air induced 
by stirring from the mixing process and then poured into 
casts. The bottom of these casts was made from glass to 
obtain smooth surfaces. The samples were cured in an 
oven at 80° C for 12 h. 

We have used sandpaper and sand-blasted PMMA as 
substrates. The sandpaper (corundum paper, grit size 
120) has the root-mean-square roughness 44 fim. From 
the measured surface topography we obtain the height 
probability distribution P(h) and the surface roughness 
power spectrum C(q) shown in Fig. [7]Ja) and[H respec- 
tively. Sand paper has much sharper and larger rough- 
ness than the counter surfaces used in normal rubber 
seal applications. However, from a theory point of view it 
should not really matter on which length scale the rough- 
ness occurs, except for "complications" such as the influ- 
ence of adhesion and fluid contamination particles (which 
tend to clog the flow channels). Nevertheless, the theory 
assumes that the average surface slope is not too large, 
and we have therefore also measured the leak rate for 
rubber seal in contact with sand-blasted Plexiglas with 
less sharp roughness. 

Our first experiment with a relative smooth Plexiglas 
(PMMA) surface showed that the leak rate decreased by 
time and finally no leaking could be observed. But this 
experiment used unfiltered tap water containing contam- 
ination particles which clogged the channels. Using dis- 
tilled water we found that the leak rate (for a given fluid 
pressure difference) to be practically time independent. 
In Fig. \T[b) and [5] we show the height probability dis- 



tribution P(h) and the power spectrum C(q) of the two 
sand-blasted PMMA used below. The root-mean-square 
roughness of the two surfaces is 34 /xm and 10 /im. 

5. Experimental results and analysis 

According to (1) and (3) we expect the leak-rate to 
increase linearly with the fluid pressure difference AP — 
P a — Pb- We first performed some experiments to test this 
prediction. In Fig. [9]we show the measured leak rate (for 
the sandpaper surface) for different fluid pressure drop 
AP for the nominal squeezing pressure Pq « 60 kPa. 
To within the accuracy of the experiment, the leak-rate 
depends linearly on AP. 

In Fig. [10] we show the logarithmic (with 10 as basis) 
of the measured leak rate for several different squeezing 
pressures (square symbols) . We show results for both the 
sandpaper surface and for the two sand-blasted PMMA 
surfaces. The solid lines are the calculated leak rate using 
the measured rubber elastic modulus E = 2.3 MPa and 
the surface power spectrum C(q) shown in Fig. O We 
show calculations using both the single-junction theory 
(blue lines) and the effective medium theory (red lines). 
In the calculations we have used a = 0.7 (for sandpa- 
per) and 0.3 (for PMMA). Note that both theories gives 
similar results for low squeezing pressures, but for larger 
squeezing pressures the effective medium theory gives a 
larger leak rate. The experimental data agree better with 
the effective medium theory than with the single-junction 
theory. 

6. Discussion 

We have presented experimental results for the leak 
rate for a PDMS rubber ring (with rectangular cross sec- 
tion), squeezed against three different surfaces: two sand- 
blasted PMMA surfaces and a sandpaper 120 surface. 
The experimental results have been compared to a simple 
single-junction theory and to a (more accurate) effective 
medium theory. The basic input in both theories is in- 
formation about the interfacial surface separation, which 
we have obtained using the contact mechanics theory of 
Persson. The pressure dependence predicted by the the- 
ory is in good agreement with the experimental data, in 
particular for the effective medium theory. 

The contact mechanics theory we use assumes ran- 
domly rough surfaces. Randomly rough surfaces have a 
Gaussian height probability distribution P(h). However, 
most surfaces of engineering interest have not Gaussian 
height probability distribution. In Fig. [7] we show the 
height distribution for the surfaces used in the present 
study. Note that P(h) is asymmetric with a tail towards 
higher h for the sandpaper surface, and towards smaller 
(negative) h for the sand-blasted PMMA surfaces. This 
is easy to understand: the sandpaper surfaces consist of 
particles with sharp edges pointing above the surface, 
while the region between particles are filled with a resin- 
binder making the valleys smoother and wider than the 
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FIG. 7: The surface height probability distribution of (a) 
sandpaper 120, and (b) two sand-blasted PMMA. The sur- 
faces have the root-mean-square roughness 44 fim, 34 fim, 
and 10 /im, respectively, and the surface area (including only 
the surface roughness with wavevector indicated in the fig- 
ure) is about 49%, 28% and 10% larger, respectively, than 
the nominal surface area Ao (i.e., the surface area projected 
on the :ry-plane). 
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FIG. 8: Surface roughness power spectrum of sandpaper 120 
and two sand-blasted PMMA. The surfaces have the root- 
mean-square roughness 44 /jm, 34 /xm, and 10 fjxci, respec- 
tively, and the surface area (including only the surface rough- 
ness with wavevector indicated in the figure) is about 49%, 
28% and 10% larger, respectively, than the nominal surface 
area Ao (i.e., the surface area projected on the xy-plane). 
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FIG. 9: Square symbols: the measured leak rate for different 
fluid pressure drop AP = P a — Pb for the nominal squeezing 
pressure Pq ~ 60 kPa. For the sandpaper surface. 
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FIG. 10: Square symbols: the measured leak rate for ten 
different squeezing pressures for sandpaper 120 (upper data 
points) and sand-blasted PMMA (lower two sets of data 
points). In each case the upper solid lines are the calculated 
leak rate using the effective medium theory and the lower solid 
lines using the single-junction theory. In the calculation we 
used the measured surface topography, the measured rubber 
elastic modulus E = 2.3 MPa and the fluid pressure differ- 
ence AP = P a — Pb = 10 kPa obtained from the height of the 
water column. In the calculations we have used a = 0.7 (for 
sandpaper) and 0.3 (for PMMA). 

peaks, which result in an asymmetric P(h) as observed. 
The PMMA surfaces are prepared by bombarding a flat 
PMMA surface with small hard particles. This result, 
at least for short time of sand-blasting, in local indenta- 
tions (where the particles hit the surface) separated by 
smoother surface regions, leading to the observed asym- 
metry in the height distribution. 

Let us now discuss how the asymmetry in the height 
distribution may effect the leak rate. To illustrate this 
we consider an extreme case: a rigid solid block with a 
flat surface in contact with a rigid substrate with periodic 
"roughness" as in Fig. QT] The substrate surfaces in (a) 
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FIG. 11: Contact between a rigid block with a flat surface 
and a rigid substrate with periodic surface structures. Two 
substrate surfaces in (a) and (b) have the same surface rough- 
ness power spectrum. Note that the empty volume between 
the surfaces is much larger in the case (a) than in case (b). 



and (b) have the same surface roughness power spectrum, 
but it is clear that in (a) the empty volume between the 
surfaces is larger than in (b), resulting in a larger leak 
rate. In the real situation the roughness is not periodic 
and the solids are not rigid, but one may expect a higher 
leak rate for the situation where the asymmetry of the 
height profile is as for the sandpaper surface. We suggest 
that this may be the physical origin of why the factor 
a is larger for the sandpaper surface as compared to the 
PMMA surfaces. Another observation which support this 
conclusion is the fact that the surface roughness power 
spectrum of the rough PMMA surface and the sandpaper 
120 surface are very similar, but the leak rate differ by 
roughly two orders of magnitude. This indicate that some 
aspects of the surface topography, not contained in the 
power spectrum, is likely to be important. We note that 
for randomly rough surface, the statistical properties of 
the surfaces are fully contained in the power spectrum 
C(q) 1 i.e., for this case only C(q) will enter in the theory 
for the leak-rate. 

To study the point discussed above, we plan to per- 
form an experiment where we "invert" the roughness of 
the sandpaper surface by producing a "negative" using 
silicon rubber. In this experiment we will squeeze a sil- 
icon rubber ring, which is cross-linked with the sandpa- 
per surface as the substrate, against a flat glass surface. 
By comparing the measured leak-rate for this configura- 
tion with that for a silicon ring with flat bottom surface 
squeezed against the same sandpaper surface, we will be 
able to address the problem illustrated in Fig. QT] 

An alternative to using the effective medium approach 
to calculate the leak-rate of seals, one may use the so 
called critical path analysis [jj. This approach has re- 



cently been applied to seals [25| but in contrast to the 
effective medium theory, there enters two parameters 
which are not easy to obtain from theory. Since the ef- 
fective medium approach has been found to be rather 
accurate (see, e.g., Ref. (26j) we believe this approach is 
more suitable for calculating the leak rate of seals. 

7. Summary and conclusion 

To summarize, we have compared experimental data 
with theory for the leak-rate of seals. The theory is based 
on percolation theory and a recently developed contact 
mechanics theory. The experiments are for silicon rub- 
ber seals in contact with sandpaper and two sand-blasted 
PMMA surfaces. The elastic properties of the rubber and 
the surface topography of the sandpaper and PMMA sur- 
faces are fully characterized. The dependence of the cal- 
culated leak-rate Q on the squeezing pressure is in good 
agreement with experiment. The simplest version of the 
theory only account for fluid flow through the percolation 
channels observed at (or close to) the percolation thresh- 
old. We have also presented another approach based on 
the effective medium approximation. This theory also 
include flow channels observed at higher magnification, 
and gives larger leak-rates than the single-junction the- 
ory, which only includes one leak-rate channel (or L y /L x 
channels for a rectangular seal). 
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Appendix A 

Here we briefly review the effective medium approach 
for calculating the fluid flow through an interface where 
the separation u(x) between the surfaces varies with the 
lateral coordinate x = (x,y). If u(x) varies slowly with 
x the Navier-Stokes equations of fluid flow reduces to 



where AP = P a — Pb is the pressure drop. 

The effective medium conductivity er c ff is obtained as 
follows. Let us study the current flow at a circular in- 
clusion (radius R) with the (constant) conductivity a lo- 
cated in an infinite conducting sheet with the (constant) 
conductivity <r e g. We introduce polar coordinates with 
the origin at the center of the circular inclusion. The 
current 

J = — cVp for r < R 

J = — (JcffVp for r > R 
We consider a steady state so that 
V- J = 

or 

V 2 p = (A4) 

If Jo is the current far from the inclusion (assumed to be 
constant) we get for r > R: 

p=[l + /(r)]J -x (A5) 

Eq. (A4) is satisfied if 

/"«+3/'(r)r- 1 =0 

A solution to this equation is / = ar~ 2 . Substituting 
this in (A5) gives 

p = [l + ar^ 2 ] J • x (A6) 

For r < R we have the solution 

P = /?Jo • x (A7) 

Since p and x • J must be continuous at r = R we get 
from (A6) and (A7): 

1 + a R- 2 = p 



J = -aVp (Al) 

where the conductivity a — u 3 (x) / '12rj . 

In the effective medium approach one replace the lo- 
cal, spatial varying, conductivity u(x) with a constant 
effective conductivity a e g. Thus the fluid flow current 
equation 

J = -CToffVp, (A2) 

as applied to a rectangular region L x x L y with the 2D 
pressure gradient dp/dx — (Pb — Pa)/L x , gives 

Q = L y J x = ^effAF (A3) 



(l — aR 2 ) cr e ff = (3a 
Combining these two equations gives 

(3 (A8) 
(Jeff + a 

The basic picture behind effective medium theories is pre- 
sented in Fig. O Thus, for a two component system, one 
assumes that the flow in the effective medium should be 
the same as the average fluid flow obtained when cir- 
cular regions of the two components are embedded in 
the effective medium. Thus, for example, the pressure 
p calculated assuming that the effective medium occur 
everywhere must equal the average c\p\ + C2P2 of the 



pressures p\ and p 2 calculated with the circular inclusion 
of the two components 1 and 2, respectively. For r < R 
we have for the effective medium p — Jo • x and using 
(A7) the equation p — c\p\ + c 2 p 2 gives 

1 = cift + c 2 (3 2 {AS) 

where c\ and c 2 are the fractions of the total area oc- 
cupied by the components 1 and 2, respectively. Using 
(A8) and (A9) gives 

2<7 e ff 2cr eff 

1 = Cl ; h C 2 ; 

Ceff + Cl CT e ff + (T2 

which is the standard Bruggeman effective medium for a 
two component system. 

If one instead have a continuous distribution of com- 
ponents (which we number by the continuous index £) 
with conductivities a = <r(£), then 

1 = y d£ P(f)/3(0 (A10) 

where is the fraction of the total surface area oc- 

cupied by the component denoted by £. The probability 
distribution P(£) is normalized so that 

J de ^(0 = 1 (ah) 

Using (A8) we get 

l=[d£P® 2 * cff (A12) 
J o- e ff + cr(ij 
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